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Abstract 

We consider a family of integer sequences generated by nonlinear re¬ 
currences of the second order, which have the curious property that the 
terms of the sequence, and integer multiples of the ratios of successive 
terms (which are also integers), appear interlaced in the continued frac¬ 
tion expansion of the sum of the reciprocals of the terms. Using the rapid 
(double exponential) growth of the terms, for each sequence it is shown 
that the sum of the reciprocals is a transcendental number. 


1 Introduction 

For some time there has been considerable interest in rational recurrences which 
surprisingly generate integer sequences. Such sequences were made popular by 
the articles of Gale Ell], who discussed some particular nonlinear recurrence 
relations of the form 


^n+N — f (^n-t-1; ■ • ■ i ^n+N— 1 ) ; (1 ■ 1) 

where / is a polynomial in A — f variables. Observe that the above recurrence 
is rational, in the sense that each new iterate Xn+N is a rational function of the 
N previous terms Xn, .. •, Xn+N-i- Starting from N initial values xq, ..., xn-i 
which are all integers, there is no reason to expect that subsequent terms will 
be, because one must divide by Xn at each step. However, a very wide variety 
of examples are now known, for which the recurrence dm has the Laurent 
property: if the initial values are viewed as variables, then for certain special 
choices of /, all of the iterates belong to the ring Z[xq ..., consisting of 

Laurent polynomials in the initial values with integer coefficients. In particular, 
the Laurent property implies that if all the initial values are taken to be 1 (or 
±1), then S Z for all n. 
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The Laurent property is a key feature of Fomin and Zelevinsky’s cluster 
algebras [5], which are generated by birational iterations of the same shape as 
cn, that is 

old variable x new variable = exchange polynomial, 

in the particular case that the exchange polynomial / is a binomial. The main 
tool available for proving the Laurent property is the Caterpillar Lemma due to 
Fomin and Zelevinsky [^, which also applies to more general choices of /, fitting 
into the broader framework of Laurent Phenomenon (LP) algebras [14] . Within 
the axiomatic setting of cluster algebras or LP algebras, there is a requirement 
that the exchange polynomials should not be divisible by any of the variables. 
However, this requirement is not necessary for the Laurent property to hold. 
Indeed, even for the case of a recurrence of second order, of the form 

^n+2 (^*^) 

the requirement that x /{f{x) is not necessary. In work by the author [lH [13] 
it was shown that recurrences of the form (lEa having the Laurent property 
ht into three classes, depending on the form of /: (i) /(O) ^ 0, in which case 
the recurrence belongs within the framework of cluster algebras (when it is 
a binomial) or LP algebras (when it is not); (ii) /(O) = 0, /'(O) ^ 0; (hi) 
/(O) = /'(O) = 0. In classes (i) and (ii) there are additional requirements on /, 
but in class (hi) one can take /(a:) = x^F{x) with arbitrary F G 1,[x\. 

The simplest non-trivial example of the form (11.21) belonging to the third of 
the classes identihed in previous work by the author |12| is the recurrence 

^n+2 (a^n+1 T !)■ (^■^) 

Due to the Laurent property, the initial values xq = xi = 1 generate an integer 
sequence: 

1,1,2,12,936,68408496,342022190843338960032,...; (1.4) 

this is sequence A112373 in the Online Encyclopedia of Integer Sequences (OEIS). 
As one might expect from the Hrst few terms, this sequence grows very rapidly: 
logXn OA" with C ~ 0.146864 and A = (3 + V^)/2. Another feature of se¬ 
quence A112373 is that the ratios yn = Xn+ifxn also form an integer sequence, 
that is 

1,2,6,78,73086,4999703411742,1710009514450915230711940280907486,..., 

(1.5) 

which is sequence A114552 in the OEIS, and the same is true for the ratios 
of ratios, i.e., = yn+i/yn = a:n-i-i -I- 1 by (11.311 : this property of the ratios 

is common to all recurrences in class (iii). Hanna made some very interesting 
empirical observations about the sequence dm 0, by considering S, the sum 
of reciprocals of the terms: 

°° 1 111 

j=0 ^ 
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In the OEIS, the digits of this number appear as sequence A114550 , yet it is not 
the decimal expansion of S that is interesting, but rather its continued fraction 
representation; with the notation 


[ao; 01,02,03,..., = oq + 


Oi 


02 


03 H- 


one finds that 

5 = [2; 1,1,2,2,6,12, 78,936,73086,68408496,4999703411742,...]. (1.7) 

What appears to be the case from the above is that (apart from the initial value 
oo = 2) the sequence of o„ in (II.7|) . which is number A114551 , is obtained by 
interlacing the original sequence (d with the ratios d- As observed by 
Shallit [16], and shown by Harris [11], this implies that the even/odd terms 
satisfy 

02n = a2n-ia2n-2> 02n+l = a2„_l(o2n + 1), (1-8) 

respectively, where the first formula holds for n > 2 and the second for n > 1. 

The purpose of this short note is to prove Hanna’s observations concerning 
the continued fraction expansion d, and generalize them to an infinite family 
of integer sequences generated by recurrences belonging to class (hi) in the 
author’s previous work [12] . At the same time we show that the number S 
given by d is transcendental, and the same is true for the sums of reciprocals 
obtained from the other sequences in this family. 

The results presented here are similar in spirit to those in a paper by Davison 
and Shallit, who found some continued fractions whose partial quotients (coef¬ 
ficients) are explicitly related to the denominators of their convergents, and 
used this to prove the transcendence of Cahen’s constant [3]. For references to 
examples of other transcendental numbers whose complete continued fraction 
expansion is known, the reader should consult the latter paper. 


2 Continued fractions for sums of reciprocals 

Before we proceed with presenting a family of sequences which generalizes d, 
we present some facts about this particular example, to motivate the proof of 
the main result. When taking the sum of reciprocals d , it is more convenient 
to exclude the index j = 0 from the sum, and then consider the partial sums 


N 


sn=j:-- 
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Calculating the finite continued fractions of these partial sums, we find that 


5l = l, = | = l + i ^3 = 1 


1 + 


1 + 


1 


1 

'+2 


5*4 = [1; 1,1,2,2,6,12], S 5 = [1; 1,1,2,2,6,12, 78,936], and 

Se = [1; 1,1,2,2,6,12, 78,936,73086,68408496] 
are the first few partial sums. As will be proved in due course, the pattern is 


Sn = N; yo,xi,yi,X2, ■ ■ ■, yN-2, a^Ar-i], (2.1) 

so that the even/odd coefficients are a 2 n = Xn and 02^+1 = 2/n respectively, and 
as we have chosen to start the sum with l/xi = 1 we now have ao = xq = 1 
which ensures that both formulae (EH hold for all n > 1. The result for S 2 
looks anomalous, but in fact (EB is seen to hold for IV = 2 upon noting that 

52 = [ 1 ; 2 ] = [ 1 ; 1 , 1 ]. 

The continued fraction for the infinite sum S^o = — is obtained in the 

limit IV —>■ 00 , and compared with (HH) we have S = Soo + 1- 

We now wish to generalize these observations to integer sequences generated 
by recurrences of the form 


Xn+2 Xn — F (^Xn+l') : 


( 2 . 2 ) 


where F{x) S Z[j:], and we assume that d = deg A > 1 to avoid a trivial case. 
If we take such a recurrence with the initial values xq = xi = 1, and set 


Vn = 


Xn+l 

Xn 


^n 


Vn+l 

y-n 


Xn+2Xn 


X 


2 

n+1 


F{Xn+l), 


then we have j/o = Ij X 2 = yi = Zo = and by induction we see that 

Xn,yn,Zn € h for all 71 > 0, SO we have three integer sequences, as long as 
the recurrence (12.2p does not reach a singularity (division by zero), which can 
happen if F{x) = 0 for some x G h. Indeed, suppose that for some positive 
integer m we have 0 G Z for 0 < n < to, but F(xm) = 0; then at the next 
step Xm+i = 0, followed by Xm +2 = 0, and then Xm+s is undefined. The analysis 
of these recurrences near to a singularity has been performed previously m- 
To avoid the possibility of reaching a singularity from the initial conditions 
xo = xi = 1, we choose F to have only positive integer coefficients {F{x) G 
Z>o[a;]), so that F{x) > 0 whenever a; > 0, and then all three sequences, (xn), 
ijjn) and (zn), consist of positive integers. In order for the continued fraction 
expansion of the sum of reciprocals to behave in the right way, we must make 
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the further assumption that F{0) = 1, which (since the degree of F is positive) 
implies i^(l) > 1, and hence Xn+i > Xn for n > 1, and the sequences (y^) and 
(zn) are strictly increasing as well. The precise rate of growth of these sequences 
will be considered in the next section, but for now we proceed with the main 
result on continued fractions. 


Theorem 2.1. For a sequence (xn) generated from the initial values xq = xi = 
1 by the recurrence h2.2\) with F(x) S Z>o[a;] and F{0) = 1, the partial sums of 
reciprocals have the continued fraction expansions 


N 


Sn = Y, 


— = [ao; ai, 02 ,..., 02^-2] 

Xj 


for all N > 1, where 


^2n — 5 


F{Xn+l) - 1 

a2n+l — - 

Xn 


G Z>o. 


(2.3) 


(2.4) 


Corollary 2.2. 

{xn) is given by 


The infinite sum of the reciprocals of the terms of the sequence 

00 ^ 

Soo — ^ ^ — [oq : Oi, 02 j ■ . ■ 5 On 5 • ■ .] j (^■^) 


i=i 


where the coefficients of the continued fraction are as in j2.4l). 


The latter result on the infinite sum follows immediately from (12.31) by taking 
the limit TV —>■ oo. Note that the convergence of the sum is guaranteed since 
the infinite continued fraction makes sense for any sequence of positive integer 
coefficients (o„); the convergence of the infinite sum can also be proved directly 
by using estimates for the growth of (a;„), as given in the next section. 


Remark 2.3. We may write 


S 


OO 


OO 




1 

ym ■■■Vj 


and observe that (y„) is a non-decreasing sequence of positive integers, with 
J/n > 2 for n> 1, so this is an example of an Engel expansion (see Theorem 2.3 
in Duverney’s book W)- 

To prove Theorem l2.ll we start by recalling some general facts about conver- 
gents of continued fractions, which are well known; for a brief summary of these 
results, the reader is referred to the first chapter of Manin and Panchishkin’s 
book m, or for more details see the book by Cassels [5]. The nth convergent 
of the continued fraction [uq; oi, 02 ,...] is given by 




where the numerators p„ and denominators are given in terms of the coeffi¬ 
cients according to the matrix identity 


ao 1 N / ai 1 \ /an 1 

1 0 J \ 1 0 J I 0 


Pn Pn-1 \ 
Qn Qn—1 J 


and both and qn are obtained recursively via the same linear three-term 
recurrence relation, that is 


Pn+l 

— O^n+lPn “t” 

Pn—1: 

(2.7) 

qn+1 

= <Xn+iqn + 

^n—li 

with the initial values 




g_i = 0, 

P-1 = go = 1, 

Po = ao. 

(2.8) 

Taking the determinant of both 

sides of (|2.6I) gives the formula 


PnQn—1 

-Pn-iqn = (- 

-1)”+', 

(2.9) 

valid for n > 0, and the identity 




PnQn—2 

-Pn-2qn = (- 

-Iran 

(2.10) 


for n > 1 follows by combining (1^ with (1^ . 

Now consider the convergents of the continued fraction whose coefficients 
a„ are given in terms of the sequence (a;„) by (12.41) . First of all, note that if 
we write F(x) = 1-1- xG(x) for G € Z>o[a:], then a 2 n+i = Xn+iG{Xn+i)/Xn = 
ynG(xn+i) G 2>o as claimed; so in general the coefficients with odd index are 
integer multiples of the ratios We prove by induction that 

q2N-l = 2/7V — 1 =-1, q2N = Xn+ 1 - (2-11) 

XN 

For N = 0 we have q-i = xi/xq — 1 = 0 and go = a:i = 1 in agreement with 
(|2.8I) . and assuming that (12.111) holds for some N, from (12.7p we have 


I / 77 / \ 1 \ XN-\-1 XN-\-1 XN-\-2 -i 

q 2 N+l = a 2 Af-|-ig 2 Ar + g 2 Af-l = (F (Xn+i) — 1)-1-1 =-1, 


XN 


XN 


XN+1 


by 


, and hence 


, fXN+2 , 

g2Ar-|-2 = a 2 N+ 2 q 2 N+l + q 2 N = Xn +1 -+ XN +1 = Xn +2 

\xn+i j 


as required. Now (12.31) is clearly true for iV = 1, and if it holds for some index 
N then, by using the second equation in (I2.11|) as well as (12.101) and (12.41) . we 
find 

o o I 1 _ P2N-2 , 1 

JN+1 = JN + 


Xn+ 1 q2N-2 


q2N 


P2Nq2N-2 — a2N + q2N-2 _ P2N 

q2N-2q2N q2N 


which means that the partial sum Sn+i is the 2iVth convergent of the continued 
fraction [cq; oi, 02 ,...] with coefficients (12.41) . so (12.31) holds for the index N+1. 
This completes the proof of Theorem 12.11 
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3 Transcendence of the sums 


We now prove the following 

Theorem 3.1. The infinite sum Sea, os in h2.5\) . is a transcendental number. 

The proof is based on Roth’s theorem, which says that if a is an irrational 
algebraic number then for an arbitrary fixed (5 > 0 there are only finitely many 
rational approximations p/q for which 


P 


a - 


q 


1 

q2+5 


(3.1) 


(see Manin and Panchishkin’s book for a brief discussion [15], and for a proof 
see the book by Cassels m)- Note that the number Sea is irrational, since its 
continued fraction expansion (ESI) consists of an infinite sequence of coefficients 
Qn ^ 0. In order to make use of Roth’s theorem, it is enough for the integer 
sequence (a:„) to satisfy the growth condition 


Xn+i > x'f for some k > 2, 


(3.2) 


for all sufficiently large n. Supposing that this is so, it follows that 

Xn+j > xf^ for j > 1 

whenever n is large enough. Then from Theorem 12.II we have 


1 P2n 
'oo 

q2n 


oo 


E 

j=n+2 


1 

- < 


OO 


E 


1 


X 


K,^ 

n+1 


Now the function g(j) = j j-i is monotone decreasing with g(2) = 2, so > jn 
for j > 2 and k > 2, which together with (I2.11|) implies 


Q P2n 
*^00 

q2n 


<E- 

1=1 


1 (1 


^n+1 


< 


^n+1 


^n+1 


K-e 

q2n 


for any e > 0 and n sufficiently large. So if e is chosen such that K — e = 2 + S>2, 
then a = Soo has infinitely many rational approximations satisfying ED), and 
hence must be transcendental. 

To show that EEj) holds for any sequence (a;„) defined by the recurrence 
(|2.2D with xq = Xi = 1 and N(0) = 1, F{x) G Z>oN as in Theorem 12.11 we 
can use a very crude estimate. Indeed, we have F{x) = 1 + ... + cx‘^ for some 
integer c > 1, so F{x) > x‘^ for all a; > 0, and then from (12.21) we obtain 

™d+2 

Xn+i > —— > ( 3 . 3 ) 

^n—1 


for all n > 1, since the ratios ?/„ = Xn+i/xn form an increasing sequence. The 
above growth condition is sufficient for p.2l) when the degree d > 2, but not 
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when d = 1, which is the case relevant to the original sequence (m. However, 

1 

this estimate can be improved upon by using Xn-i < in the first inequality 
in (13.311 ■ to yield 

X„+l>Xn^"~^ >x^„/^ (3.4) 

for d > 1 . 

In fact, we can get a more accurate measure of growth from asymptotic 
arguments. Upon taking the logarithm of ( 1 ^ we find that A„ = log Xn satisfies 

A„+i - (d + 2)An + A„_i = logc + a„, with = log ( , ( 3 , 5 ) 

V cxi J 

Note that a„ = log(l + 0{x~^)) = as n —>■ 00 , which means that, to 

leading order, the growth of A„ is governed by a homogeneous linear equation 
with constant coefficients, given by the vanishing of the left-hand side of (13.511 . 
The characteristic equation is — (d -I- 2)A -1-1 = 0, with the largest root being 


^ ^ d + 2 + d(d -I- 4) ^ 2 


(3.6) 


Thus we find that 


A„ ^ CA", 


(3.7) 


for some C > 0. Hence A„+i/A„ —)> A as n —>■ 00 , and so for any e > 0 it follows 
that A„+i > (A — e)A„ for all sufficiently large n, giving the required growth 
condition 

^n+l ^ ^n 

The asymptotic properties of the sequence (a;„) can be determined more 
precisely by adapting the methods of Aho and Sloane [I] , leading to the following 
result, which is easy to verify directly from (13.511 . 

Proposition 3.2. For the initial conditions xq = xi = 1, the logarithm An = 
logXn, of each term of the sequence satisfying 112. is given by the formula 


(l-A-i)A"-(l-A)A- 


A„, = ( -1 j logcd-t^ (:____— j ak, (3.8) 




^n — k _ ^k—n 


where ak is defined as in (TOP and A as in Eg). 


Corollary 3.3. To leading order, the asymptotic approximation of the logarithm 
An is given by fd. 7|), where 


\-\ 

and for the terms of the sequence 




-k 


O^k, 


k^l 


c-^exp(CA^). 









Remark 3.4. The form of the expression is the discrete analogue of the so¬ 
lution of a linear inhomogeneous differential equation as obtained via the method 
of variation of parameters m; if ak were given in advance, then it would pro¬ 
vide the exact solution of ^3.5\] , viewed as a linear equation with initial values 
Ao = Ai = 0. However, as was pointed out for a different example by Aho and 
Sloane m, a formula such as iTOI) only represents the solution of the corre¬ 
sponding nonlinear equation (in this case the equation h2.0i) for Xn) in a tau- 
tologous sense, because ak depends explicitly on the terms of the sequence (x„). 
Nevertheless, this formula does yield useful asymptotic information about the 
sequence. 
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